HRI-P-12-01-002 



Path Integral Junctions 



Satoshi Ohya 

Harish- Chandra Research Institute 
Chhatnag Road, Jhusi, Allahabad 211 019, India 

E-mail : ohyaOhr i . res . in 

(Dated: June 6, 2012) 

Abstract 

We propose path integral description for quantum mechanical systems on compact 
graphs consisting of N segments of the same length. Provided the bulk Hamiltonian is 
segment-independent, scale-invariant boundary conditions given by self-adjoint exten- 
sion of a Hamiltonian operator turn out to be in one-to-one correspondence with N x N 
matrix- valued weight factors on the path integral side. We show that these weight fac- 
tors are given by ./V-dimensional unitary representations of the infinite dihedral group. 



1 Introduction 



It has been long known that, in path integral on multiply-connected space A4, linear com- 
bination coefficients (weight factors) for homotopically distinct sectors are given by unitary 
representations of the fundamental group iri(M). For standard "scalar" quantum me- 
chanics, i.e. the case where wave functions are given by single-component scalar-valued 
functions, these are one-dimensional unitary representations [IHH and hence characterized 
by U(l) phase. As first discussed by Horvathy et al. [5j, however, for "vector" quantum 
mechanics, i.e. the case where wave functions become multi-component vector-valued func- 
tions, weight factors could be higher-dimensional unitary representations [5| and hence in 
principle enjoy much more parameters than those of "scalar" quantum mechanics case. 

In the present paper we attempt to formulate path integral description of boundary 
conditions provided by Kirchhoff's law of probability current [6] (self-adjoint extension of 
Hamiltonian operator) on quantum graphs, by continuing our previous investigation on the 
path integral on star graphs [TJ. The key to understanding the path integral description of 
self-adjoint extension parameters is higher-dimensional unitary representations. Originally, 
Horvathy et al. focused on internal symmetries as the origin of vector- valued wave functions 
[5]. However, the origin of vector- valuedness is not necessarily internal symmetries. An 
alternative realization of vector- valued wave functions is folding trick [U|9], which maps a 
defect system to a boundary system. In this paper we will work on a compact graph that 
consists of N segments of the same length £, (0,1) 3 Xi (i = 1, • • , N), each of which is 
connected at xi = and £, respectively; see figure [TJ Such graph can be folded into a 
system on a single segment with A-component vector-valued wave functions, with U(N) 
family of boundary conditions at each boundary. An important point to note is that this 
graph contains a large number of topologically distinct compact graphs as its "subgroups" ; 

namely, a reduction U(N) ->■ U(M\) x • • • x U(M n ) (Mi H YM n = N) leads to a number 

of compact graphs that have different topologies, which is nothing but the topology change 
discussed in fT0j[TTJ . Typical examples of such compact graphs are depicted in figure [2j 

A puzzle may arise from homotopy structures in folded theories, however. Typically, 
folded systems boil down to systems on half- line or segment, which are simply-connected 
spaces such that its fundamental groups are trivial. Nonetheless, path integral on, say, a 
segment is known to be given by a sum of partial amplitudes corresponding to an infinite 
number of classical bouncing paths (for textbook exposition see e.g. chapter 6 of |12j). 
Actually, this is an evidence that such system also intrinsically possesses its own group 
structure in the path space. 
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Figure 1: Junctions of N segments of the same length i. Dashed lines represent the current 
flow among different segments. Arrow indicates the direction of coordinates. 
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Figure 2: Typical segment junctions and its foldings, (a) the case of N = 2, (b) the case of 
the reduction U{N) t -> [U(l) e ) N , and (c) the case of the reduction U(2N) e -> [U(2) e ] N . 



The purpose of the paper is twofold. Firstly, we would like to clarify the group structure 
of the space of all possible paths on a segment. Although the path integral on a segment has 
been studied over the years [T3HT6] . to the best of our knowledge, little is known about its 
group structure. We show that the path space on a segment is furnished with the structure 
of the infinite dihedral group = 7L x Z2. Secondly, we discuss that path integral on 
a segment junction is given by a linear combination of partial amplitudes, each of which 
corresponds to a distinct sector. We show that, for the case of scale-invariant subfamily 
of boundary conditions, the weight factor for each sector falls into an iV-dimensional 
unitary representation of D^, which is the main result of the present paper. 

The organization of the paper is as follows: We will begin in section [2] by studying a sys- 
tem for a free spinless particle on segment junctions in the folding picture, which is exactly 
solvable and exemplifies enough the path structure of the system. This non-interacting the- 
ory provides the simplest yet concrete illustration for the group theoretical structure of the 
weight factors. In section [3] we generalize to interacting theories by introducing segment- 
independent bulk interactions. We show that the weight factors in the path integral are 
generally given by iV-dimensional unitary representations of the infinite dihedral group . 
Section H] is devoted to conclusions. Computational details are relegated to appendix [Al 

Throughout the paper we will work in the units h = 2m = 1. 
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2 Scale-invariant boundary conditions 



Let us first start with quantum mechanics for a free particle on junctions of N segments 
of the same length I. By folding trick [U[9], such system is mapped to an alternative 
equivalent system on a single segment (0,£), whose Hilbert space is given by a tensor 
product H = L 2 (0,£) <g> C*. There, the wave function becomes an TV-component vector- 
valued function ip(x) = (ipi(x), ■ ■ ■ ,^n(x)) t G H; see figure [TJ In this system self-adjoint 
extension of free Hamiltonian IHf ree = — -j^z <8> In leads to the following U(N)q x U(N)i 
family of boundary conditions |17j : 

(I N - Uo)^(0) - iL(I N + Uo)^(0) = 0, U G U(N) , (2.1a) 
(I N - U £ )tp(£) + iL(I N + V t )$(£) = 0, U, G U(N) t , (2.1b) 

where In is an N x N identity matrix, L is an arbitrary length scale and ip'(x) = ^-(x). The 



sign difference ±iL is just for later convenience. We note that figure 2(a), 2(b), and 2(c) 
correspond, respectively, to the cases of N = 2, U(N) e [U(l) e ] N , and U(N) t -> [U(2)^ 1 
(N: even integer). 

We now focus on the case of scale-invariant boundary conditions, which are specified by 
hermitian unitary matrices 0[18] that satisfy Uq = U| = In- in that case the boundary 
conditions (|2.1aj) (|2. lbj) are simply reduced to the following conditions: 

(Ijv - Uo)^(0) = (I N + Uo)^(0) = 0, (2.2a) 
(I N - V^(£) = (I N + U e )$(£) = 0, (2.2b) 

which follow from multiplying equations (|2.1ap and (|2. lb>|) by Ijv iUo and In ±Uf, respec- 
tively. Note that the boundary conditions (|2.2ap and (j2.2b|) do not involve the scale parame- 
ter L, as it should. Now it is easy to solve the Schrodinger equation Wf TCC ip(x) = Eip(x) and 
find the complete orthonormal set of energy eigenfunctions with respect to (|2.2ap and (|2.2bj) . 
By expanding the kernel in terms of the energy eigenfunctions and then using the argument 
principle, which enables us to switch from the summation over energy spectrum to the sum- 
mation over (classical) trajectories, the total Feynman kernel ~K(x,y;T) = (x|e _iHflccT |y) is 
cast into the following N x N matrix- valued form (see appendix : 

K(x, y; T) = ^(U U,)^ free (x, y - 2n£; T) 

+ ^(V V e ) n V K iTee (x,-y-2n£-T), 0<x,y<£, (2.3) 

whose ij-component gives the transition amplitude for a particle propagating from the 
position y on the jth segment to the position x on the ith segment. Kf Tee (x, y; T) is the scalar 
Feynman kernel for a single free particle on a whole line R and given by Kf ree (x,y;T) = 
-j= = exp [ij (^y^) ]. Notice that the scalar kernel is invariant under the action of both 
reflection and continuous translation, Kf ree (— x, —y\ T) = Kf Tee (x, y; T) and Kf Tee (x + a,y + 
a;T) = Kf ree (x,y;T) (a: arbitrary constant). We also note that the total Feynman kernel 
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(|2.3|) satisfies the following conditions: 



(Initial condition) K(ac, y;0) = 5(x — y)Ijv; (2.4a) 

(Composition rule) / dz K(x, z; Ti)K(z, y; T2) = K(x,y;Ti +T2); (2.4b) 

Jo 

(Unitarity) K?(x, y; T) = K(y, x; — T); (2.4c) 

(Boundary conditions) (I N - U )K(0, y; T) = (I N + V )(d x K)(0, y; T) = 0, (2.4d) 

(1 N - V e )K(£, y; T) = [l N + V e ) (d x K) (£, y; T) = 0. (2.4e) 



We now turn to the physical interpretation of the total Feynman kernel (|2.3p . The 
exponent ^ ^ x-(±y-2ni) j - g j^i^g j-, u t ac tion integral 5[x c i] = J^dtji^t) for a free 
particle on R traveling along the classical trajectory x c i(t) = X{ + (^jr^-)t with initial and 
final positions Xj = ±y — 2n£ and Xf = x. Thanks to the translation- and reflection- 
invariance of the scalar kernel, these straight line trajectories on R can be interpreted as 
bouncing paths on the segment (0,^); see figure^ These bouncing paths are weighted by 
the scale-invariant S-matrix Uo (U^) every time when a particle bounces off the boundary 
x = (x = £). (Uo and play the roles of boundary S-matrices; see appendix 1X1 ) 
For example, for a classical straight line trajectory with initial position y — 2n£ (n > 0), 
which corresponds to the bouncing path that hits the boundaries 2n times in alternating 
order, the weight factor is given by an alternating product of hermitian unitary S-matrices 



UoVe ■ ■ - VoVe = (UoU^) ra ; see figure 3(a) for the case of n = 1. Likewise, for a classical 
trajectory with initial position — y — 2n£ (n > 0), which corresponds to the bouncing path 
that hits the boundaries 2n + 1 times in alternating order, the weight factor is given by 
UoU^ • • -UoU^Uo = (VoVe) n Vo; see figure [3(b)] for the case of n = 1. The summation over 
all integer n thus ensures the summation over all possible bouncing paths on the segment. 

3 Weight factor and Af-dimensional unitary representation 

of 9* Z 2 * Z 2 

So far we have seen that the hermitian unitary matrices Uo and U^ which specify the scale- 
invariant boundary conditions on the operator formalism side have appeared as weight 
factors on the path integral side. In this section we first clarify the group structure of the 
path space on a segment and then show that, without invoking the boundary conditions, 
Uo and U^ are introduced via A^-dimensional unitary representations of the infinite dihedral 
group. This algebraic consideration leads to the generalization to systems in the presence 
of bulk interactions. 

To begin with, let us first recall basics of the infinite dihedral group = Z x Z2, which 
is generated by two operations, a discrete translation T € Z and a reflection 1Z G Z2, whose 
actions on R can be defined as T : y 1— > y — 2£ and 1Z : y 1— > —y. All elements of are then 
of the forms T n and T n lZ(= lZT~ n ) with integer n, which act on R as T n : y 1— > y — 2n£ 
and T n lZ : y 1— > —y — 2n£, and obey the following multiplication rules: 

r n or m = r n+m , T n oT m n = T n+m n, (s.ia) 

V'TZ o T m = T n ~ m n, T n KoT m K = T n ~ m , (3.1b) 
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(c) Quantum fluctuation of (a) (weight factor: U^UfUoUf = UoUf 
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(d) Quantum fluctuation of (b) (weight factor: UoU^U^U^Uo = UoUfUo) 



Figure 3: Typical classical trajectories and its quantum fluctuations. 



where circle (o) stands for abstract group multiplication. Note that identity element is T°, 
and inverse elements are (T n ) -1 = T~ n and (T n TZ)~ 1 = T n TZ. Note also that T n in ^ 0) 
and T m TZ are non-commutative, T n o T m TZ 7^ T m TZ o T n ; that is, Z?oo is a non-Abelian 
discrete group. 

Let us next clarify the group structure of the space of all possible paths on a segment. 
As we have seen in section [21 the total Feynman kernel on a segment is given by a linear 
combination of partial amplitudes. Each partial amplitude is obtained by integration over 
one particular path class, which is a set of paths on M that cross each wall between the 
initial and final positions odd number of times. Such classification of paths is furnished with 
the infinite dihedral group structure. Indeed, as a natural generalization of the classical free 
particle trajectories, paths shall acquire the S-matrix Un (U^) every time when a particle 
crosses the walls at x = —2n£ (x = — (2n + 1)^); see figure 3(c) and |3(d)[ Any paths 
shall therefore be labeled by certain product of S-matrices. Then, thanks to the hermitian 
unitarity of the S-matrices, Uq = Ijv and \jf =In, any weight factor falls into an alternating 
product of Uo and U^. Thus, without any loss of generality, we may denote these path classes 
as [(UoU^)"] and [(UoU£) n Uo]. Composition of these path classes are given by, for example, 
[U U,U ] * [U U,] = [UnU.UnUoUf] = [U ] and [U U,] * [U U,U ] = [U U,U U £ U ] . By 
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making use of the fact that is isomorphic to the free product Z 2 * Z2, which we will 
revisit shortly, we can easily check that these path classes {[(UoU^)" - ], [(UoU£) n Uo]} satisfy 
the -Dqo = 7* 2 * Z2 multiplication rules (|3.1a|) (|3.1bp under the composition of paths. A 
crucial point to note is that these path classes are simply classified by its initial positions 
T n y = y — 2n£ and T n 1Zy = —y — 2n£. 

Now we are ready to discuss the weight factor in path integral on segment junctions. 
As the simplest generalization of free particle case discussed in the previous section, we 
will concentrate on a system described by the segment-independent bulk Hamiltonian that 
takes the form Hbuik = -Htmik <8> Ijv. In this case a kernel that acts on % = L 2 (0,£) (g) C N 
should be factorized as (scalar kernel) x (A x N constant matrix). As discussed in the 
previous paragraph, the total Feynman kernel must be given by a linear combination of 
partial amplitudes, each of which is weighted by an A x N matrix- valued linear combination 
coefficient. Putting these considerations together, the total Feynman kernel must be of the 
form 

K(x, y;T) = Y J [W(T n )K R (x, T n y; T) + W(T*Tl)K*(x, TKy\ T)], 0<x,y<£, (3.2) 

where {W(7 _ri ), W(T n TZ)} are N x A constant matrices. K^(x,y;T) is the scalar kernel on 
M. and assumed to satisfy the following conditions: 

(Initial condition) K^(x,y;0) = 5(x — y); (3.3a) 

/oo 

dz K K (x, z; Ti)K R (z, y; T 2 ) = Kr(x, y; Ti + T 2 ); (3.3b) 
-00 

(Unitarity) K&(x, y; T) = K R (y, x; — T); (3.3c) 

(Doo symmetry) K R (x,y;T) = K^{Qx,Qy;T), \/Q G D^. (3.3d) 

We note that fl m symmetry will be justified by translation- and reflection- invariant bulk 
scalar Hamiltonian, T'.ffbulkT" = ^-^buik^ -1 = -f^buik- Below we will show that, under 
the assumptions (|3.3ap - (|3.3dp . the ansatz (|3.2p will satisfy the initial condition (|2.4a|> . com- 
position rule (|2~4R and unitarity O^cj) if and only if the weight factors {W(T n ), W(T n K)} 
are given by A-dimensional unitary representations of -Doo- After proving this statement, 
we will rederive the scale-invariant boundary conditions by using unitary representations. 

The proof is almost parallel to that presented in the previous work [7] such that we here 
discuss only the composition rule. By substituting the ansatz (|3.2p . the left hand side of 
(|2.4bp becomes 

2 W(T n )W(T m ) I dzK w {x,T n z-,T l )Ki S .(z 1 T m y-,T 2 ) 

+ J] W(T n Tl)W(7 m n) I dzKK{x,T n Uz-T l )Ki ll {z,T n lly-T 2 ) 

+ V W(T n )M{T m K) I dzK*(x,T n ^T 1 )K*(z,T n nv\T 2 ) 
^» Jo 

+ V W(T n U)W{T m ) (dzK w (x,T n nz-T x )K^z,T m y^T 2 ). (3.4) 
_tt^ J 
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Noting that the identities K R (z,T m y;T 2 ) = K R (T n z,T n+m y;T 2 ) and K R (z,T m Ky;T 2 ) = 
K^(T n TZz, T n ~ m y; T 2 ), which follow from symmetry (|3.3d|) and multiplication rules 
dSJajl plbjl . the first two integrals in fl33} become f* dz K R (x, T n z; Ti)K R (z, T m y; T 2 ) = 
f^^d^K^x^^K^T^y^) and f£ dz K^(x,T n TZz; Ti)K^(z, T m 1Zy; T 2 ) = 

J-{2n+V\l d£ C; 7i)^r(C T n ~ m y; T 2 ), where we have changed the integration variables 

as £ = T n z = z — 2n£ and £ = T n TZz = —z — 2n£, respectively. Similarly, the last 
two integrals in (J33J can be written as f^ +%)e d£K m (x, £; Ti)K*(£, T n+m Ky; T 2 ) and 
J_(2n+i)^ d£,K^(x, £; Ti)ifR(^, T n ~ m 1Zy; T 2 ). Collecting the above pieces the left hand side 
of (I2.4bp is cast into the following form: 



E 



/•(-2n+l)(. t—2ni 

W(T n )W(T m ) / dt + WiT^WiT^K) / d£ 

J-2n£ J-(2n+l)e 



•(-2n+l)£ j—2ni 

-(2n+l)£ 

X K v (x,(;T 1 )K x (S,T n+m y,T 2 ) 



+ E 



n,mdi 



W(T n )W(T m K) / d£ + W(T n K)W(T- m ) / df 



x K M (x,^,T 1 )K^,r n+m 'Jly,T 2 ), (3.5) 

where in the second and fourth terms we have changed the summation variable as m — > —m. 
Notice that the integral /l 2 n™ +1 ^ ^ anc ^ f-(2n+i)i ^ cover om y the shaded and unshaded 
regions of figure [31 respectively. Thus, in order to realize the right hand side of (|2.4bj) with 
the assumption (j3.3bj) . W must satisfy the following conditions: 

W(T n )W(T m ) = W(T n K)W(T- m K) = W(T n+m ), (3.6a) 
W(T n )W(T m K) = W(T n K)W{T- m ) = W(T n+m K), (3.6b) 

which are nothing but multiplication rules ()3.1ap and (|3.1bp . Hence W must be an Nx N 
matrix representation of D^. Similarly, one can obtain the conditions for the weight factors 
by imposing the initial condition (|2.4ap and the unitarity (|2.4c[) . It is straightforward to 
show that these two conditions lead to the following constraints for the weight factors: 

W(T°) = In, (3.7) 

and 

W f (T n ) = W(T~ n ), W\T n K) = W(T n K), (3.8) 

from which one deduces that the weight factors must be unitary matrices, W _1 (T n ) = 
W+CT 1 ) and W- 1 (T n K) = W^(T n TZ). Putting all these results together, we find that W 
must be an TV-dimensional unitary representation of the infinite dihedral group 

W : Doo -»■ U(N). (3.9) 

Now, the problem is reduced to the problem of how to construct such unitary represen- 
tations. To solve this, we first note that the infinite dihedral group is isomorphic to the free 
product of two cyclic groups of order 2, Z 2 *Z 2 , each of which is generated by two distinct 
parity transformations Vq and Vt that satisfy Vq = Vf = 1. Defining the actions of these 
two parity transformations on R as Vq '■ y >• — y and Vt '■ y i— > 2£ — y, whose fixed points 
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are y = and y = £, respectively, we find that Vq and VqPi correspond to the generators 
1Z and T. Noting that every element of Z2 * Z2 is given by an alternating product of Vo 
and Vi, we immediately see that the map Z2 * Z2 — > given by {PoPtf 1 1— > T n and 
(T~ , o'Pe) n 'Po T n TZ. is an isomorphism. The second point to note is that iV-dimensional 
unitary representation of Z2 is just given by an N x N hermitian unitary matrix [7j. There- 
fore, assigning two distinct hermitian unitary matrices to Vq and Ve, we obtain the following 
iV-dimensional unitary representations of D^: 

W(T") = (U iy n , (3.10a) 
W(T n TZ) = (U iy n U , (3.10b) 

where Uo,U^ G U(N) are hermitian unitary matrices that satisfy Uq = = In- These 
weight factors are nothing but those obtained in the case of free particle. 

Now it is straightforward to derive the boundary conditions for the total Feynman kernel 
K. If the scalar kernel is continuous and smooth at x = and £, which is the case of 
iCfree, we have 



K(0,y;T) = [In + W(1Z)] W(T n )K$i(0, T n y; T), (3.11a) 

(d x K)(0,y;T) = [I N - W(K)} ^ W(r n )(a^ R )(0, T n y; T), (3.11b) 

K(£, y; T) = [I N + W(ftT)] ]T W(T n )K R (£, T n y; T), (3.11c) 

(d x K)(£,y;T) = [In — W(TZT)] ^ W(T n )(d x K]&)(£, T n y; T), (3.11d) 

which follow from W(RT n ) = W{K)W(T n ), W{TZT n+1 ) = W{KT)W(T n ) and sym- 
metry (|3.3dp . Making use of the unitary representations W(1Z) = Uo and W(TZT) = Ug, 
we immediately reproduce the scale-invariant boundary conditions (In — Uo)IK(0, y; T) = 
(I N +Vo)(d x K)(0,y;T) = and (J N -W e )K(l,y;T) = (I N +V e ){d x K)(£, y; T) = 0, where we 
have used the orthogonality relations of hermitian unitary matrices (In ± Uo)(Iat =F Uo) = 
(I JV ±U^)(I 7V T^) = 0. 

Path integral representation of the total Feynman kernel (|3.2p is now easy. Since path 
integral representation of is a standard problem, we just assume that is given by 
the standard configuration space path integral. Then, if the bulk interactions are segment- 
independent, the path integral for systems on junctions of iV segments of equal length £ is 
given by the following N x N matrix-valued functional integral: 

rx(T)=x / ,-T \ 

K(x,y;T)= V W(<?) / Vx(t)exp (i / dt L(x(t),x(t)) ) , (3.12) 
geDoo Jx(o)=g y V Jo J 

where W(Q) (Q G -Doo) is given in (l3~TUaD (l3~T0bl . L(x(t),x(t)) is a generic .Doo-invariant 
one-particle Lagrangian, and Dx(t) is a -invariant path integral measure on R. 

4 Conclusions 

In this paper we have studied path integral description of quantum mechanics on junctions 
of N segments of equal length £. We have seen that scale-invariant subfamily of boundary 
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conditions is well translated into the problem of weight factors in path integral. Provided 
the bulk Hamiltonian is segment-independent, the weight factors are generally given by 
iV-dimensional unitary representations of the infinite dihedral group = Z2 * Z2; that 
is, in this case the problem of boundary conditions boils down to the problem of unitary 
representations. These two problems are in one-to-one correspondence with each other, 
since both two are solved once two hermitian unitary matrices Uo and U^ are specified. 

It should be noted that our result ()3. 12j) does not cover the whole self-adjoint domain 
of Hamiltonian operator. For non-scale-invariant case, the boundary S-matrices generally 
depend on the particle momentum such that the total Feynman kernel cannot be factorized 
as (scalar kernel) x (TV x N constant matrix) in the configuration space path integral. As 
discussed in the previous work [7] , the key to understanding the non-scale-invariant bound- 
ary conditions would be to generalize the weight factors in the phase space path integral. 
We hope to address this issue elsewhere. 



Appendix A Feynman kernel for a free particle 

In this appendix we first solve the Schrodinger equation with scale-invariant boundary 
conditions (|2.2aj) (|2.2bp and then derive the Feynman kernel (|2.3p by using the complete 
orthonormal basis of energy eigenfunctions and the argument principle. 

The general solution to the Schrodinger equation Wf rce ip(x) = Eip(x) for positive energy 
E > is given by a linear combination of plain waves 

J{x; p) = A e~ ipx + A e e- ip{i - x \ p = Ve > 0, (A.l) 

where Aq and are iV-component coefficient column vectors. Substituting the general 
solution to the boundary conditions (|2.2a|) (|2,2b|) we get the following eigenvalue equation: 

MM g) - g) , (A,) 
where M(p) is a 2iV x 2iV unitary matrix given by 

M W = U o) e '" (A ' 3) 

It is worth mentioning here the physical meaning of the hermitian unitary matrices Uo and 
U^. By using the relations A^ = Uoe ip Vlo and Aq = Uge ip Ag, which follow from f|A.2|) . 
the solution (jA.lj) can be rewritten as ${x\p) = (I N e~ ipx + U e ipx )A = {l N a ipi - x ~^ + 
U ee -ip(x-t))A £ . These expressions imply that Uo and U^ play the roles of boundary S- 
matrices at x = and £, whose r/-component gives the transmission amplitude (or reflection 
amplitude if i = j) for a particle traveling from the jth segment to the ith segment; see 
figure |H 

Let us next discuss the complete orthonormal set of energy eigenfunctions. The con- 
dition ()A.2p implies that the eigenvectors and eigenvalues of M(p) will respectively deter- 
mine the coefficient vectors Aq, Ai and the energy spectrum via quantization condition 

(eigenvalue of M(p)) = 1. Let e lf)a (0 < 9 a < 2tt) and A a = (5 ) (a = 1, • • • , 27V) be the 

ath eigenvalue and its corresponding normalized eigenvector of the unitary matrix ( u ^ ) 
that satisfy the eigenvalue equation v e )A a = e tea A a , orthonormality A^A^ = 5 ab and 
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Imp 



-X X X X- 



IPI Pi Pi ••■ 

-* X X X- 



Rep 



x = 



x = t 



Figure 4: One-particle boundary S-matrices 
U and U £ . 



Figure 5: Simple poles (x) and integration 
contour C. 



completeness X*a=i A a A aJf = I2N- The normalized energy eigenfunction is then given by 
$ a (x;p%) = ^(lge-^« x + A^e-^-^) with p a n = > (n = 1, 2, • • • ) being the 

positive roots of the quantization condition e l P e -+ iea = 10 These normalized energy eigen- 
functions provide the complete orthonormal basis of the system such that the Feynman 
kernel is expanded as K(x,y;T) = E™=i ^ fl (3;;K) e ~ !T ' P "' (Vi Pn) j which can be 

written as follows: 



2N oo 

K(x,y;T) = - £ J2E 

a=l n=l 



e ip « x I N 
e ipl{t-x)l N 



N 



(A.4) 



where P a = A a A a ^ is a hermitian projector onto the ath eigenspace for the unitary matrix 
M(p). (Note that the energy eigenfunction can be written as ip a (x;p) = ^{jpy-xfj ) A a .) 

Now we wish to switch from the summation over the energy spectrum to the summation 
over (classical) trajectories, which can be achieved via the following contour integral (i.e. 
argument principle): 



oo „ 
7^1 J c 



^/(P)^log(l-e^)) 
C 2m dp 



E 

m=l 



m=0 



£^0 + 



0+ie 

oo 



0-ie 



-im(p£+O a ) 



! W" ^ e^+^/fr), (/: test function), 



(A.5) 



where the integration contour C is chosen to enclose all the positive roots of the quantization 
condition in a counterclockwise direction; see figure [5j The first equality just follows from 
the residue theorem and the second equality the geometric series expansion: 



d 



i(p£+6 a ) 



m=l 

oo 

m=0 



,im(p<+fl») for Imp >0, 
r im(p«+fl") for imp <0 . 



Since there are no longer poles in the integrand, we can take the limit e 
line of ([A~5l) . 



(A.6) 



0_i_ in the third 



1 For the sake of simplicity we here assume that 6 a ^ for any a, which guarantees the absence of constant 
zero-modes in the spectrum. For the case of the presence of zero-modes, one has to modify the integration 
contour C in figure[5]by changing the leftmost vertical line to an infinitesimal semicircle. However, the result 
(I A, 71) remains unchanged regardless of the presence or absence of zero-modes. 
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Now, with the aid of the formula (|A.5p the Feynman kernel (|A.4|) is cast into the following 
form: 



lilt- 



dp ( e^Ijy 



t 



-iTp 2 



(A.7) 



which follows from the spectral decomposition M. m (p) = Y^l 1 e Lrn ( pi+Ba ^ a . By using the 
definition (|A.3|) . the mth power of unitary matrix M(p) is easily calculated with the result 



M m (p) = < 



OW 1 \ &ip2ne 
v (UoU^)™^ 

o (u,u )"uA eJp(2n+1 

(U U,)"U j 



for m = 2n, 



for m = 2n + 1 , 



(Ai 



where n is an integer. Substituting (|A.8j) into (|A.7[) one immediately sees that the total 
Feynman kernel is recast into the following form: 



K(x, y; T) = £(U U*) n f°° ^ exp fiT 



x — y + 2n£ \ > 
P I m ~P 



tip 

,2^ 



exp ( iT 



x + y + 2ni i .-, 
p 



(A.9) 



By performing the momentum integration, we obtain the result (|2,3[) . 
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